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INTRODUCTION 



Short-range-interaction spin-glass (SG) models have raised a lot of contro- 

versies in the last decades. From the theoretical point of view, the case of Ising spins 
represents the most convenient to be investigated, in such a way that a large effort has 
been dedicated to the understanding of the Ising SG model. The majority of works 
concentrated on three-dimensional Ising SG models, for which, besides its physical real- 
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ly accepted nowadays that a SG phase occurs at finite temperatures 

3,Q- nnnn 

hierarchical lattices have been a very useful tool for the study of SG models [5|, la, [ij, llSl . 
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26|, essentially due to the possibility of carrying out 



exact calculations, or performing relative^ 



Apart from a few exceptions 
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low-time-consuming numerical computations. 
261], most of the hierarchical lattices considered 
so far belong to the Migdal-Kadanoff (MK) family. Taking into account the significant 
reduction of efforts in the investigation of these models, some of the results obtained are 
quite impressive: (i) The lower critical dimension, below which the SG phase transition 
occurs at zero temperature, which was the object of a lot of controversy about three 
decades ago, was correctly estimated within a MK renormalization-group (RG) approach 
[sl, almost a decade before the consensus that this quantity should be greater than 2, 
)ut smaller than 3, through different numerical approaches, like numerical simulations 
and zero-temperature domain- wall arguments [8|]; (ii) The critical-temperature 
estimates of Ref. [sl for an Ising SG on a MK hierarchical lattice of fractal dimension 
D = 3, with symmetric distributions, are [kBTd J) = 1.05 ± 0.02 (±J distribution) and 
{kBTc/ J) = 0.88 ± 0.02 (Gaussian distribution of width J). The most recent Monte 
Carlo simulations on a cubic lattice [l^ yield {ksTc/J) = 1.120 ± 0.004 in the first 
case, whereas for the later, {ksTc/ J) = 0.951 ± 0.009, leading to relative discrepancies 
of about 4% when compared with the results of Ref. [5|, taking into account the error 
bars; (iii) A zero-temperature analysis of a Gaussian Ising SG, on a special hierarchical 
lattice with fractal dimension D = 2 yielded an estimate for the stiffness exponent 
y {y = —l/u, where u is the exponent associated with the divergence of the correlation 
length at zero temperature) in agreement with those obtained from other, more time- 
consuming, numerical approaches; (iv) The same hierachical lattice of Ref. [l7|, mentioned 
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above, produced a precise ferromagnetic-paramagnetic critical frontier for the ±J Ising 
SG model ^ . 



A major question in the SG problem nowadays concerns the applicability of some 
results from the mean-field solution for short-range-interaction systems. In particular, 
whether the SG phase is properly described by an infinite number of order parameters 
(i.e., an order-parameter function {27]), manifesting the property of replica-symmetry 
breaking (RSB) [l^; moreover, if an Ising SG, in the presence of an external magnetic 
field, exhibits the Almeida-Thouless (AT) line 28j, which separates a low-temperature 
region characterized by RSB, from a high-temperature one, described in terms of a single 
order parameter, along which the replica-symmetric solution holds. Numerical simula- 
tions for nearest-neighbor-interaction Ising SGs are always hard to be performed, since 
equilibration becomes difficult for large system sizes and low temperatures. However, in 
spite of the small lattice sizes considered, there are evidences from Monte Carlo simula- 
ions that a critical line in the presence of a field exists in d = 4 29|, but not in d = 3 
111 . [l3 |. Furthermore, a zero-temperature analysis of the energy landscape in the case 
(i = 3 is compatible with a transition from the SG to the paramagnetic phase for a finite 
critical field, although the possibility of a critical field equal to zero was not excluded 30 1. 
However, it is possible to have a critical frontier separating the SG and paramagnetic 
phases, for low-dimensional short-range-interaction Ising SGs in the presence of a mag- 
netic field, which is not an AT-like line. In order to ensure that this critical line represents 
a true AT line, one should also verify evidences of RSB below such a frontier. This possi- 
bility would correspond to an "intermediate" scenario [sO], between the mean-field RSB 
solution and the much simpler droplet picture 0.0.1 In this case, one would expect 
that RSB effects should appear below this line, at some finite dimension d, contrary to 
claims of the droplet model, according to which the AT line should occur only in the limit 
of infinite dimensions. 

One of the advantages for the study of SGs on hierarchical lattices is that one does not 
go through equilibration difficulties. In view of this, the D = 3 MK hierarchical lattice has 
been used also for an investi gati on of RSB in the low-temperature phase on an Ising SG 



without a magnetic field 18, 
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|2J]. For temperatures in the range 0.7 Tc < T < Tc 



a picture showing characteristics of RSB was observed, although for lower temperatures 
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FIG. 1: Basic cells of three hierarchical lattices that approach the cubic lattice, (a) The MK cell 
of fractal dimension D = 3 (usually called of diamond cell); (b) The dual of the diamond cell 
(fractal dimension D = 3); (c) The tridimensional Wheatstone-Bridgc cell (fractal dimension 
D ~ 3.58). The empty circles represent the external sites, whereas the black circles are sites to 
be decimated in the renormalization process. 

the results agree with the simpler, replica-symmetric scenario. 

To our knowledge, short-range-interaction SGs in the presence of an external magnetic 
field have never been investigated on hierarchical lattices. In the present work, we study 
an Ising SG model, in the presence of different types of external magnetic fields, on 
three hierarchical lattices (two of them with a fractal dimension D = 3 and another one 
characterized hj D ^ 3.58). In the next section we define the model, the hierarchical 
lattices, and the numerical procedure. In section III we present and discuss our results. 

II. THE MODEL AND NUMERICAL PROCEDURE 

Herein, we study an Ising SG, in the presence of an external magnetic field, 

n^-Y,•J^JS^S,-Y,h,S, {S, = ±l), (1) 
<ij> i 

where {Jij} denote random couplings between two spins located at nearest-neighboring 
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sites i and j of a given hierarchical lattice, following a symmetric Gaussian probability 
distribution, 

P(J,,) = -L= exp (-^] . (2) 



v/2^ " V 2J2 
For the magnetic fields we consider three different cases, namely, 

P{h,) = 6{hi - Ho) (uniform field), (3) 



P{hi) = - 5{hi — Hq) + - 5{hi + Hq) (symmetric bimodal distribution), (4) 



Pihi) = r- — ^ exp 



2a2 



^2 



(Gaussian distribution). (5) 



The Hamiltonian of Eq. ([T]) will be investigated on three different hierarchical lattices 
that approach the cubic lattice. These lattices are generated by starting the process from 
the 0th level of the lattice-generation hierarchy, with a single bond joining the external 
sites (denoted by fi and u). Then, in each iteration step one replaces a single bond by 
a unit cell, like the ones shown in Fig. [T], in such a way that in its first hierarchy, each 
lattice is represented by a unit cell; the hierarchical lattice is constructed up to a given 
Af-th hierarchy {Af ^ 1). The cells in Figs. [D^a) and (b) are considered as dual to one 
another, and their results for critical temperatures usually represent lower and upper 



limits, respectively, with respect to the correct value of the cubic lattice [3l(|. The cell 
of Fig. Wic) is a three-dimensional Wheatstone Bridge cell, and to our knowledge, it has 
never been used in the study or random magnetic models. 

The RG procedure works in the inverse way of the lattice generation, i.e., through a 
decimation of the internal sites of a given cell, leading to renormalized quantities associ- 
ated with the external sites. Defining the dimensionless couplings and fields, Kij = (3Jij , 
Hi = phi [P =JJ{kBT)], the corresponding RG equations may be written in the general 



form (see Ref. 



32| for the explicit form of these equations for the MK cell of Fig. [U^a)), 



6 



1 . (Z^^Z_ 



+ 



4 ° V 4 



where Zs^^s^ represent partition functions of a given cell with the external spins kept fixed 

%,5. = Tr|5^ (i^^,^)} [exp(-/?7^)] . (9) 

It should be noticed that the Hamiltonian associated with any of the cells of Fig. [T] may 
be split into H = H' + H^S^ + H^jS^, where 7i' represents the Hamiltonian of the cell 
with Hfj^ = Hu = 0. Therefore, the partition function Zs^^Si, niay be rewritten as, 

Zs^,s. = eMH.S, + H,S,)Zs^,s^ , (10) 



where 



2s^,s. = Tr{s, ii^f.,u)} [expi-pn')] (11) 



represents the partition function of a given cell without considering the contributions of 
the fields on its external sites. Accordingly, Eqs. ([S])^(IH]) take the following form. 



<. = 7log(^^^) , (12) 



< = ^. + ilog(|±±|±^) , (13) 
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Thus, we clearly see through Eqs. (fT3|) and (fT^ that after the RG transformation, the 
field on each remaining site depends on the previous field on this site plus a contribution 
due to the decimation of the inner spins in the cell. 

If at the beginning of the RG procedure (A/'-th hierarchy), one uses the probabil- 
ity distribution of Eq. ([2]) for the couplings, and one of the probability distributions 
of Eqs. ([Sj-Q for the fields, one has at this level, the average values, < Kij >= 0, 
< Hi >= Hq (uniform field and Gaussian distribution), or < Hi >= (symmetric bi- 
modal distribution). An important quantity to be used herein is the ratio of associated 
widths, 

r = ^ ; aK=< {K,^- < K,, >)' ■ an =< {H~ < H^ >f >^l^ , (15) 

which, at the AA-th level, is infinite for a uniform field (since ou = in this case), 
whereas r = JjH^ (bimodal distribution) and r = J/a (Gaussian distribution for the 
fields). However, as the renormalization procedure goes on (in fact, right after the first 
RG transformation), the system of coupled equations that define the renormalization 
[Eqs. dSD^dS])] introduces correlations between the couplings and fields, and as a result of 
this, one has a joint probability distribution, PjKij, Hi, Hj), to be followed. 

We have used the method proposed in Ref. [3^1] to follow this distribution numerically. 
This technique consists in generating a pool of M triplets {Kij, Hi, Hj}, initially chosen 
by generating random numbers according to the above distribution for the couplings, and 
one of the probability distributions for the fields. An iteration consists in M operations, 
where in each of them one picks randomly triplets from the pool (each chosen triplet 
is assigned to a bond in one of the cells of Fig. [1]) in order to generate the effective 
quantities of Eqs. ([6])-([8]) that will define a triplet of the renormalized pool. At zero 



temperature Eqs. ([6])-([8]) become much simpler (see, e.g., Ref. 3^ for the explicit form 



of these equations for the MK cell of Fig. [11(a)), and we have used these simplified forms 
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FIG. 2: Evolution of the ratio r with the RG steps n, for an Ising SG on the hierarchical lattice 
defined by the unit cell of Fig. ^a), in the presence of an initial uniform field Hq. (a) Typical 
values of Hq/J at zero temperature; the critical field is in the range 0.52 < (Hqc/J) < 0.54. 
(b) Typical temperatures for {Hq/J) = 0.3; the critical temperature is in the range 0.65 < 
{kBTjJ) < 0.70. 

to determine the zero-temperature critical points. It is important to notice that in this 
procedure there may occur superpositions of random fields in given sites of the unit cells; 
whenever this happens, we consider an arithmetic average of the superposed random fields. 
After each iteration, one calculates the lowest moments associated with the couplings and 
fields, and in particular, the ratio r defined in Eq. f|T5l) . 

Below, we present and discuss the results obtained from this formalism. 

III. RESULTS AND DISCUSSION 

For the results that follow, we have used M = 160000, although we checked that our 
results did not change (within the error bars) for larger pools of triplets. For all three 
hierarchical lattices considered, the first moments presented small fluctuations around 
their initial values [i.e., < Kij >= and < Hi > (either zero or not)] under successive 
renormalizations. Therefore, our critieria for the identification of the attractors, and their 
associated phases, was based on the widths ax, cth, and their ratio r, which showed that 
the behavior of ctx always prevailed over the one of an- 

We have found typically two distinct behaviors for these quantities, under successive 
RG transformations, which we associated with two distinct phases, as described below. 
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FIG. 3: Phase diagrams for an Ising SG on the hierarchical lattice defined by the unit cell of 
Fig. [U^a), in the presence of an initial random magnetic field following a Gaussian probability 
distribution, (a) Plane Hq versus temperature, for different values of a. (b) Plane a versus 
temperature, for Hq = 0. P and SG denote the paramagnetic and spin-glass phases, respectively. 

(i) Paramagnetic (P) phase: this occurs for sufficiently large temperatures and/or fields, 
where ax ^ 0, an ^ 0, with r ^ 0. (ii) SG phase: this phase appears for low 
temperatures and fields, being characterized by ax oo and r oo. In most of the 
cases, we have found throughout this phase an increase on an as well, but still keeping 
r — > oo; however, in some situations (essentially for symmetric field distributions), we 
have found that ax — ^ oo and (Xh — * 0. The critical frontier separating these two phases 
was considered as the one where the parameter r changes very slowly. Typical behaviors 
of the ratio r, under successive RG transformations, are illustrated in Fig. [2] for an Ising 
SG on the hierarchical lattice defined by the unit cell of Fig. [T](a), in the presence of an 
initial uniform field Hq. In Fig. [2t^a) we present this quantity for different values of Hq/J 
at zero temperature, whereas in Fig. [2](b), we have fixed the initial field [{Hq/J) = 0.3] 
and have varied the temperature; in both figures one sees that the ratio of widths r may 
present qualitatively distinct behaviors, depending on the initial parameters considered, 
signalling different attractors of the renormalization. 

In Fig. [3] we exhibit phase diagrams obtained from the present RG approach for the 
Ising SG defined by Eqs. ([I]), ([2]), and ([5]), on the diamond hierarchical lattice of Fig. [T](a). 
In Fig. [3]^a) we present phase diagrams in the plane Hq versus temperature for different 
choices of a; one notices that for increasing values of a, the SG phase decreases; a similar 
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FIG. 4: Phase diagram for an Ising SG on the hierarchical lattice defined by the unit cell of 
Fig. dja), in the presence of an initial random magnetic field following a symmetric bimodal 
probability distribution. 

effect has been observed in mean-field theory, where this phase is delimited by an AT 



line, associated with RSB [S^]- It is important to stress that the case a = in Fig. [3]^a) 
corresponds, within the RG procedure, to an initial uniform field; therefore, the uniform 
and Gaussian random fields produce a qualitatively similar critical frontier. In Fig.[3](b) we 
show the phase diagram in the plane a versus temperature, for the case of a symmetric 
Gaussian distribution for the fields {Hq = 0), which displays also a SG phase for low 
temperatures. 

In Fig. m we exhibit the phase diagram for the present Ising SG model under the 
bimodal random field defined by Eq. (jlj). One notices that the critical frontier of Fig. H] 
presents a change of concavity for low temperatures and so, it is qualitatively different 
from those of a uniform field [case cr = in Fig. [3]^a)] and Gaussian random field [cases 
(T > in Fig.[3t^a)]. Curiously, this critical frontier is very similar to the one of a symmetric 
Gaussian distribution, shown in Fig.[3t^b). Analogous to the previous cases, shown in Fig. [3] 
(uniform and Gaussian random fields), the most important aspect, i.e., the existence of 
the SG attractor, associated with a SG phase at low temperatures, applies for the bimodal 
random field, as well. 
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We have also investigated the present Ising SG model on the hierarchical lattices of 
Figs-lU^b) anddl^c); although the results are quantitatively different, the qualitative behav- 
iors of the phase diagrams are the same as those shown above for the diamond hierarchical 
lattice. In Tables I and II we compare critical parameters, associated with the phase dia- 
grams, for the three hierarchical lattices investigated. In Table I we present the values of 
the critical temperatures for Hq = 0, in the case of the uniform and bimodal distributions 
for the fields, whereas for the Gaussian distribution, the results refer to the particular 
width (cr/J) = 1. In Table II we present values of critical fields, at zero temperature; for 
a Gaussian distribution, our zero-temperature results correspond to the cases, (cr/J) = 1, 
or Hq = (cf., e.g.. Fig. [3]). /^From these tables one notices that the critical parameters 
for the hierarchical lattices defined by the cells in Figs. [l](a) and (b) are always below 
and above, respectively, the ones of the hierarchical lattice of Fig. [I^c). This confirms the 
current belief that these lattices yield lower and upper limits for phase-diagram critical 
parameters [3l|; in addition to that, it suggests that the Wheatstone-Bridge lattice of 
Fig. [T](c) may lead to good approximations in the study of three-dimensional SG systems. 
In fact, the corresponding critical-temperature estimate for Hq = (uniform and bimodal 
cases) shown in Table I, {ksTc/J) = 0.980(1), is in good agreement with the recent esti- 
mate from Monte Carlo simulations on a cubic lattice, {ksTc/ J) = 0.951(9) [l^, leading 
to a relative discrepancy of 2%, taking into account the error bars. 

The results above are reinforced by a zero-temperature analysis of the evolution of the 
widths associated with the couplings (uj) and fields {cfff), with respect to the RG steps 



aj ~ ; an ~ 5™ , (16) 



which are characterized by the exponents y (usually known as the stiffness exponent) 
and u] the scaling factor is 6 = 2 for all three cells of Fig. [H In Fig. [5] we illustrate the 
power-law behaviors of Eq. (fT6l) for the case of an initial uniform field Hq on a hierarchical 
lattice defined by the unit cell of Fig. [TJ^a). The values of Hq chosen in Fig. [5] correspond to 
zero-temperature points in the spin-glass phase of the case a = in Fig. [3t^a). The scaling 
range (in the RG steps n) starts for small n (typically around n = 5) for < ifo ^ 
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FIG. 5: Evolution of the widths associated with the couplings (crj) (empty symbols) and fields 
{(th) (black symbols) with the RG steps n for typical values of an initial uniform field Hq. In 
the first case we exhibit also the case Hq = 0. 

but for larger values of Hq (essentially as one approaches the critical value Hqc) the scaling 
forms of Eq. ( |T6l) are only satisfied for higher ranges of the iteration steps n. We have 
noticed that, in all cases, the exponents y and u are weakly dependent on the initial value 
Hq (as can be seen on the examples shown in Fig. [5]); as a consequence of this, we have 
assumed universality with respect to Hq, for each of these exponents. Our estimates, for 
the case of an initial uniform field Hq, are presented in Table III for each of the hierarchical 
lattices of Figs. [TJ the error bars take into account slight variations in these exponents 
throughout the interval < Hq < Hqc. The fact that y > u supports the existence of a 
low-temperature SG phase, based on the criterion described above for the quantity r of 
Eq. (USD. 

At this point, it is important to stress that for the families of hierarchical lattices 
associated with the cells in Figs. [TJ^a) and (b), for which the fractal dimension may be 
changed easily by varying the number of parallel paths, we have also investigated lattices 
with fractal dimensions D < 3 (more specifically, those characterized by three parallel 
paths, for which D ^ 2.585). In these cases, we did not find any evidence of a SG attractor, 
suggesting that for these lattices one has a lower critical dimension, 2.585 < Di < 3. 
Taking into account the above-mentioned properties of these lattices, concerning lower 
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Hierarchical 


ksTc/J 


ksTc/J 


lattice 


Uniform and Bimodal 


Gaussian [(cr/J) = 1] 


Cell 1(a) 


0.880(1) 


0.661(1) 


Cell 1(b) 


1.761(1) 


1.750(1) 


Cell 1(c) 


0.980(1) 


0.891(1) 



TABLE I: The critical temperatures for Hq = 0, with the fields following Eqs. ([3])~(l5]), for the 
hierarchical lattices defined by the cells of Fig. [TJ In the case of the Gaussian distribution for 
the fields, we have chosen (cr/J) = 1. The error bars refer to the usual approach to criticality 
characteristic of the RG technique. 



Hierarchical 


Uniform 


Bimodal 


Gaussian [(cr/J) 


= 1] 


Gaussian {Hq = 0) 


lattice 


Hqc/ J 


Hqc/ J 


Hoc/J 




CTc/J 


Cell 1(a) 


0.530(1) 


3.919(1) 


0.433(1) 




3.036(1) 


Cell 1(b) 


1.414(2) 


21.45(2) 


1.418(1) 




20.25(2) 


Ceh 1(c) 


0.590(2) 


5.907(2) 


0.566(2) 




5.884(2) 



TABLE II: Zero-temperature critical values of the fields following Eqs. ([3])-([5]), for the hierar- 
chical lattices defined by the cells of Fig. [TJ In the case of the Gaussian distribution for the 
fields, we have chosen either (cr/J) = 1, or Hq = 0. The error bars refer to the usual approach 
to criticality characteristic of the RG technique. 

and upper limits for the critical temperatures, one may expect that these bounds for 
the lower critical dimension, associated with a SG phase in the presence of an external 
magnetic field, should apply to other lattices as well. 

To conclude, we have investigated a nearest-neighbor-interaction Ising spin glass model, 
in the presence of an external magnetic field, on three different hierarchical lattices that 
approach the cubic lattice. In the beginning of the renormalization-group procedure, 
the magnetic field was considered as uniform, or randomly distributed, following either 
a bimodal or a Gaussian probability distribution. In all cases considered, a spin-glass 
attractor was found, in the plane magnetic field versus temperature, which was associated 
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Hierarchical 


Exponent y 


Exponent u 


lattice 






Cell 1(a) 


0.245(8) 


0.140(8) 


Cell 1(b) 


0.254(2) 


0.002(1) 


Cell 1(c) 


0.223(2) 


0.007(3) 



TABLE III: The zero-temperature exponents y and u (defined in Eq. (jl6p ) for an initial uniform 
field Hq in the range < Hq < Hqc, where Hqc is given in Table II for each of the hierarchical 
lattices defined by the cells of Fig. [TJ 



with a low-temperature spin-glass phase. In the particular cases of hierarchical lattices for 
which the fractal dimension may be changed easily by varying the number of parallel paths, 
we have verified that the lower critical dimension, associated with a finite-temperature 
spin-glass phase, lies in the interval 2.585 < Di < 3. The present results show that, in 
what concerns the hierarchical lattices studied, there is a spin-glass phase in the presence 
of an external magnetic field, contrary to the claim of recent numerical simulations on 
the cubic lattice [ill, Since the above results concern the existence of a spin-glass 
attractor, under renormalization-group transformations, the question of replica-symmetry 
breaking throughout this phase, which would characterize the critical lines presented 
herein as Almeida- Thouless lines, represents a point that deserves further investigations. 
Taking into account that hierarchical lattices have been a useful tool for studying spin- 
glass systems, the present results motivate the investigation of replica-symmetry breaking 
properties for spin glasses on hierarchical lattices, in the presence of external magnetic 
fields. It is possible that the picture found pre vioulsy on the Z) = 3 Migdal-Kadanoff 
hierarchical lattice, without a magnetic field [l8|, |2l|, |22, l24] , exibiting replica-symmetry 
breaking characteristics only close to the critical temperature, may change under the 
presence of an external magnetic field. 



Acknowledgments 



15 



We thank Prof. E. M. F. Curado for fruitful conversations. The partial financial supports 
from CNPq and Pronex/MCT/FAPERJ (Brazilian agencies) are acknowledged. 



[1] V. Dotsenko, Introduction to the Replica Theory of Disordered Statistical Systems (Cam- 
bridge University Press, Cambridge, 2001). 

[2] H. Nishimori, Statistical Physics of Spin Glasses and Information Processing (Oxford Uni- 
versity Press, Oxford 2001). 

[3] Spin Glasses and Random Fields, edited by A. P. Young (World Scientific, Singapore, 1998). 

[4] K. H. Fischer and J. A. Hertz, Spin Glasses (Cambridge University Press, Cambridge, 
1991). 

[5] B. W. Southern and A. P. Young, J. Phys. C 10, 2179 (1977). 

[6] R. N. Bhatt and A. P. Young, Phys. Rev. Lett. 54 , 924 (1985). 

[7] A. T. Ogielski and I. Morgenstern, Phys. Rev. Lett. 54 , 928 (1985). 

[8] A. J. Bray and M. A. Moore, in Heidelberg Colloquium on Glassy Dynamics, Lecture Notes 
in Physics 275, Eds. J.L. van Hemmcn and I. Morgenstern (Springer- Verlag, Berhn, 1987). 

[9] R. N. Bhatt and A. P. Young, Phys. Rev. B 37, 5606 (1988). 
[10] E. Marinari, G. Parisi, and J.J. Ruiz-Lorenzo, in Spin Glasses and Random Fields, edited 

by A. P. Young (World Scientific, Singapore, 1998). 
[11] A. P. Young and H. G. Katzgraber, Phys. Rev. Lett. 93, 207203 (2004). 
[12] H. G. Katzgraber, M. Korner, and A. P. Young, Phys. Rev. B 73, 224432 (2006). 
[13] T. Jorg, H. G. Katzgraber, and F. Krz^kala, Phys. Rev. Lett. 100, 197202 (2008). 
[14] J. R. Banavar and A. J. Bray, Phys. Rev. B 35, 8888 (1987). 
[15] J. R. Banavar and A. J. Bray, Phys. Rev. B 38, 2564 (1988). 
[16] M. Nifle and H. J. Hilhorst, Phys. Rev. Lett. 68, 2992 (1992). 
[17] F. D. Nobre, Phys. Lett. A 250, 163 (1998). 

[18] M. A. Moore, H. Bokil, and B. Drossel, Phys. Rev. Lett. 81, 4252 (1998). 

[19] E. Nogueira Jr., S. Coutinho, F. D. Nobre, and E. M. F. Curado, Physica A 271, 125 



16 



(1999). 

[20] E. M. F. Curado, F. D. Nobre, and S. Coutinho, Phys. Rev. E 60, 3761 (1999). 

[21] E. Marinari, G. Parisi, J. J. Ruiz-Lorenzo, and F. Zuliani, Phys. Rev. Lett. 82, 5176 (1999). 
[22] H. Bokil, A. J. Bray, B. Drossel, and M. A. Moore, Phys. Rev. Lett. 82, 5177 (1999). 
[23] F. D. Nobre, Physica A 280, 456 (2000). 

[24] B. Drossel, H. Bokil, M. A. Moore, and A. J. Bray, Eur. Phys. J. B 13, 369 (2000). 

[25] F. D. Nobre, Phys. Rev. E 64, 046108 (2001). 

[26] F. D. Nobre, Physica A 319, 362 (2003). 

[27] G. Parisi, Phys. Rev. Lett. 43 , 1754 (1979); 50, 1946 (1983). 

[28] J. R. L. de Almeida and D. J. Thouless, J. Phys. A 11, 983 (1978). 

[29] D. Badoni, J. C. Ciria, G. Parisi, F. Ritort, J. Pech and J. J. Ruiz-Lorenzo, Europhys. Lett. 
21, 495 (1993). 

[30] F. Krzakala, J. Houdayer, E. Marinari, O. C. Martin, and G. Parisi, Phys. Rev. Lett. 87, 

197204 (2001). 
[31] J. R. Melrose, J. Phys. A 16, 1041 (1983). 
[32] A. Rosas and S. Coutinho, Physica A 335, 115 (2004). 
[33] M. S. Cao and J. Machta, Phys. Rev. B 48, 3177 (1993). 

[34] R. F. Soares, F. D. Nobre and J. R. L. de Almeida, Phys. Rev. B 50 , 6151 (1994). 



